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Abstract 

We present a complexity reduction algorithm for a family of parameter-dependent linear systems when the 
system parameters belong to a compact semi-algebraic set. This algorithm potentially describes the underlying 
dynamical system with fewer parameters or state variables. To do so, it minimizes the distance (i.e., iioo-norni 
" of the difference) between the original system and its reduced version. We present a sub-optimal solution to 

this problem using sum-of-squares optimization methods. We present the results for both continuous-time and 
discrete-time systems. Lastly, we illustrate the applicability of our proposed algorithm on numerical examples. 

9-i 1 Introduction 

1 Large-scale systems are often composed of several interacting subsystems described by local parameters that need 
C**) to be identified when designing model-based control laws. The parameters are typically a function of the working 
' points of the subsystems and their physical properties. Hence, they vary over time based on the operation mode. 

In practice, we like to develop a family of controllers that only depend on a few of the system parameters, such 
\ that we do not need to adjust the whole controller whenever a parameter changes in the system. In addition, we 
jyy . might want to study the relative importance of the system parameters. Gain scheduling and supervisory control 
^ | are examples of parameter-dependent controllers (THZ] - However, these design methods implicitly assume in most 
O ■ cases that the overall controller has access to the entire set of model parameters. This assumption might not be 
realistic in many practical cases (see [5] and references therein for a detailed discussion). Hence, we are interested in 
introducing a complexity reduction algorithm to effectively remove some of the system parameters or to decrease its 
^ ' order while preserving the input-output transfer function to some extent. Doing so, we can then simplify the control 
r**~ ■ design procedure or satisfy the requirements described above on model parameter dependencies. The problem of 
^T*" " model reduction for parameter-dependent linear systems has been studied extensively [9TU3j. For instance, the 
authors in |10) used a multidimensional system formulation and introduced a generalization of controllability and 
observability Gramians using a pair of linear matrix inequalities (LMIs). Using these generalized Gramians, they 
performed balanced truncation to extract the reduced system. They also calculated an upper bound for the error of 



o 
in 



Ov 

1 this truncation. However, the reduced order system presented in [TU| is not optimal since the introduced upper bound 
for the truncation error is not necessarily tight. In this paper, we introduce a near-optimal numerical procedure for 
extracting these reduced system. 

Specifically, we are interested in minimizing the iJoo-norm of the difference of the original system transfer functions 
1 and its reduced version over a compact semi-algebraic set of system parameters. Using the bounded real lemma |14j . 
^ . we transform this problem to a parameter-dependent feasibility-checking bilinear matrix inequality (BMI). We use 
the method of alternating LMIs [T5] to transform this parameter-dependent BMI into a string of LMIs. Then, we 
use the method introduced in |16j to solve these parameter-dependent LMIs by means of sum-of-squares optimiza- 
tion. This algorithm results in a sub-optimal solution because (1) when using the method of alternating LMIs, we 
cannot guarantee the convergence of the proposed algorithm (i.e., there exists always a BMI such that you cannot 
check its feasibility using the method of alternating LMIs [15]), and (2) when using sum-of-squares optimization for 
solving the parameter-dependent LMIs, the lack of convergence to a solution does not imply the infeasibility of the 
original problem (since a sum-of-squares matrix is indeed a positive-definite polynomial matrix but not the other 
way around) [16]. Due to relying on sum-of-square optimization, the proposed algorithm does not scale well with 
the system dimension and the number of parameters. However, we might be able to exploit sparsity patterns or 
symmetry structures in future to develop better numerical algorithms |17j . Despite these inefficiencies, we observe 
that the proposed algorithm is fairly strong in solving the proposed numerical examples in Section 2] 

Recently, there have been many studies on using sum-of-squares optimization methods in control design |18H21| . 
For instance, the problem of finding a polynomial Lyapunov function for nonlinear systems was considered in [181119] . 
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The problem of optimal linear quadratic control design for parameter-dependent discrete-time systems was discussed 
in |20j . However, to the best of our knowledge, no attention has been paid to complexity reduction using sum-of-square 
optimization. 

The rest of the paper is organized as follows. In Section[2J we present the mathematical problem formulation. We 
introduce our complexity reduction algorithm and prove its suboptimality in Scction[3J Wc illustrate the applicability 
of the proposed algorithm on two numerical examples and compare their results with available methods in Sectional 
Finally, we conclude the paper in Section [51 

1.1 Notation 

The sets of integer, natural, real, and complex numbers are denoted respectively by Z, N, K, and C. For any 
integer number n G Z and any real number x G M, wc define the notations Z>(>)„ = {m G Z | m > (>)n} and 
IR > (>) :r — {y G K | y > {>)x}, respectively. All other sets are denoted by calligraphic letters such as A and B. 

Matrices are denoted by capital roman letters such as A and B. A > (>)0 means that the symmetric matrix 
A G R" xn is positive definite (positive semidefinite) and A > (>)B implies that A — B > (>)0. 

The ring of polynomials with coefficients in R is denoted by K[a], where a is the vector of variables. For any 
given n,m G N, a polynomial matrix X(a) G M.[a] nxm is a matrix whose entries are polynomials in R[a], that is, 
Xij(a) G R[a] for 1 < i < n and 1 < j < m. For any given n G N, a matrix polynomial X{a) G R[a] nx ™ is positive 
definite (positive semidefinite) if for each a £ A, the matrix X(a) G R Ttxn is positive definite (positive semidefinite), 
where the set A will be defined in the text. 

For any given n G N, a matrix polynomial X(a) G R[a]" xn is a sum-of-square matrix, denoted by X(a) >z 0, 
if there exits a matrix polynomial Y(a) G R[a] nx ™ such that X(a) — Y(a) T Y(a). We introduce the notation 
S[a] n — {X(a) G M[a] nx ™ | X(a) >z 0} to capture the set of all sum-of-square matrices. When n = 1, we use S[a] 
instead of S [a] 1 . 



2 Problem Formulation 

In this section, wc present the mathematical formulation of the complexity reduction problem introduced in Section [1] 
for both continuous-time and discrete-time parameter-dependent linear systems. 

2.1 Continuous-Time Systems 

Consider a parameter-dependent continuous-time linear dynamical system described by 

r(*-„\ ■ jx{t) = A{a)x{t)+B{a)u(t), m 
1 ' ' ' \y(t) - C(a)x(t) + D(a)u(t), [L) 

where x(t) G M™ is the state vector, u(t) G K m is the control input, y(t) G M° is the system output, and a G W is 
the parameter vector. Note that in ([1]), we use the notation 

G(s; a) = C(a)(sl - A{a)y 1 B{a) + D(a). 

Throughout this paper, wc assume that a G A C M p , where A is defined to be the set of eligible parameters. We 
are interested in extracting a reduced parameter-dependent continuous-time linear system described by 

r'(«- n >\ ■ j*'(t) = A'(a')x'(t) + B'(a')u(t), 

1 ' ' ' \y'{t) - C'(a')x'(t) + D'(a')u(t), [Z) 

where x'(t) G M™ is the reduced system state vector, y'(t) G R° is its output, and a' G A' C W is the reduced 
parameter vector. Note that the output vector dimension stays the same. We define the reduced set of eligible 
parameters as 



A' = \ a' G R p ' I 3£ G 



a' 



eA 



Remark 2.1 We name this procedure as complexity reduction because we can potentially reduce the number of the 
parameters with which the system is described (since, by definition, we assume p' < p). In addition, by choosing 
n' < n, we may also reduce the system order. 

Throughout this paper, we make the following assumption concerning the model matrices: 
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Assumption 2.2 The model matrices in (QP and are polynomials in terms of the system parameters a and a' , that 
is, A{a) G R[a] nxn , B{a) G R[a} nxm , C{a) G M[a] ox ™, D(a) G R[a] oxm , A'(a') G R[a!'] n ' xn ' } B'(a') G K[a'] n ' xm , 
G'(a') G R[a'] oxn \ and D'(a') G R[a'] oxm . 

We are interested in finding G'(s; a') to minimize the distance between the systems in ([I) and ([2]): 

inf sup \\G(s; a) -G'&Tia))^, (3) 

G'(s;q') „ e _4 

where the projection T : W — > MP is defined as T(x) = [x\ ■ ■ ■ x p /] T for all iel p . The optimization problem in 
is to be solved subject to the reduced system state-space description (|2|) and the fact that the model matrices are 
polynomial matrices in a and a' (Assumption |2?2|) . 



Remark 2.3 If we are ultimately interested in designing a controller using the reduced system, we should solve the 
optimization problem 

inf sup \\G(s;a)- 1 (G(s;a)-G'(s;T(a)))\\ oo , 

see J22j/ . In the case that G(s; a) does not vary much over the set of eligible parameter A, we can instead solve the 
optimization problem 

inf sup || G(s; (G(s; a) - G'(s;T(a))) \\ (4) 

for some fixed ft G A.. As our developed algorithm would not change much for solving instead of we would 
only focus on solving (0) in this paper. 

2.2 Discrete-Time Systems 

Consider a parameter-dependent discrete-time linear time-invariant system 

ru- n \ ■ jx(k + l) = A(a)x(k) + B(a)u(k), 
Uy ' a> ■ \ y{k) = C(a)x(k) + D(a)u(k), 

where, similar to the previous subsection, x(k) G M n is the state vector, u(k) G K m is the control input, y(k) G M° is 
the system output, and a G M p is the parameter vector. In we use the notation 

G( 3 ; a) = C(o)(aJ - A(a)y 1 B(a) + D(a). 

We define the reduced system as 

C'U- n'\ • \ X '^ k + J ) = A '(*'>'( k ) + B'{a')u(k), 

u ^ a >-\ y'(k) = C'(a')x'(k) + D'(a')u(k), [> 

where x'(t) G M™ is the reduced system state vector, y'(t) G R° is its output, and a' G A' C K p is the reduced 
parameter vector. For these parameter-dependent discrete-time systems, we are interested in solving the optimization 
problem 

inf sup ||G( 3 ; a)- G'^^a))!^ , (7) 

G'{y,a') a ^A 

subject to the reduced system state-space description in @ and Assumption 12.21 In the next section, we present 
solutions to the optimization problems ^ and ©. 

3 Main Results 

In this section, we rewrite the optimization problems as parameter-dependent feasibility-checking BMIs. We use the 
method of alternating LMIs to transform this parameter-dependent BMI into a string of LMIs, which we then solve 
using sum-of-squares optimization methods. First, we present the solution for continuous-time systems. 
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Procedure 1 Extracting the sub-optimal reduced system G'(s;a'). 



Input: 5 G R> , e G R> , 7 G R> , n' G N, p' G N, p' G N, d P G N, eN, d B e N, d c eN,d D e N, and d Qf G N 
for all g < I < L. 

Output: A'(a') G R[a'] n ' x ™', P'(a') G R[a']"' xm , C'(a') G R[a'] ox "', D'(a') G R[a'] 0Xm , and P(a) G 5[a]™+™'. 
Initialization: Pick A'(a') G R[a'] ,l ' x, \ B'(a') G R[a']"' xm , C'(a') G R[a'] oX "', and D'(a') G M[a'] oxm such that 
deg(A'(a')) = d A , deg(B'(a')) = d B , deg(C"(a')) = d c, and deg(D'(a')) = do, respectively. Also, set P id(a) = 0. 



for f = 1,2,... do 

Find polynomial matrix P(a) G <S[a] n+ " with deg(P(a)) = dp and polynomial matrices Qe(a) G S[a] n+n + m +° 
with deg(Q£(a)) = dQ e for all 1 < £ < L, such that 



A(a) T P(a) + P(a)A(a) * * 
B{a) T P(a) -I * 

(5(a) 5(a) -7 2 / . 



(P.i) 



Find polynomial matrices Qi(a) G S[a] n+n + m +° with deg(Q^(a)) = dQ e for all 1 < £ < L and model ma- 
trices A'(a') G R[a']"' xn ', B'(a') G R[a'] Tl ' xm , C'(a') G R[a'] ox "', and D'(a') G R[a'] oxm with respectively 
deg(A'(a')) = d A , dcg(P'(a')) = d B , deg(C*'(a')) = d c , and deg(D'(a')) = d D such that 



A(a) T P(a) + P(a)A(a) * * 
B(a) T P(a) -I * 

(5(a) 5(a) -7 2 / 



+ Qo + ^Q£(a)g«(a) = 0. 
( -1 



(P.2) 



if max Qe ^ \\P(a) - P id(a)|| < $ then 

break 
else 

Poid(a) <- P(a) 
end if 
end for 



3.1 Complexity Reduction for Continuous-Time Systems 

Before stating the results, let us define the augmented system as 



x(t) 
±'{t) 



A(a) 

y(t)-y'(t) = C(a) 



x(t) 
x'(t) 

x(t) 

At) 



B(a)u(t), 
D(a)u(t). 



where 



and 



A(a) 



A{a) 
A' {a') 



, B{a) 



B(a) 
B>(a>) 



C(a) = [ C(a) -C'(a') 1 , D(a) = D(a) - D'{a'). 



(8) 



(9) 



(10) 



Now, we are ready to present the first result of the paper. The next lemma transforms the iJoo-optimization problem 
in ^ into a parameter-dependent BMI. 

Lemma 3.1 For a fixed a G A and G'(s;a'), we have \\G(s;a) — G'(s; c/)||oo < 7 */ and only if there exists 
P{a) = P(a) T G R(«+«')x («+«') such that p( a ) > an d 



A{a) T P{a) + P{a)A{a) * * 
B(a) T P(a) -I * 

C(a) 5(a) -7 2 / 



< 0, 



for all a G A. 



Proof: The proof follows from Bounded Real Lemma [T3] on the augmented system 
a G A and G'(s;a'), the augmented system is simply a linear time-invariant system. 



(11) 



Note that after fixing 
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Note that Lemma 13.11 does not guarantee that P(a) is a matrix polynomial in a. In the next lemma, we show 
that this is indeed the case using the results in [23] . 

Lemma 3.2 Let A be a compact subset of M 11 . Then, for a fixed G' (s; a'), ||C?(s; a) — G'(s; a/JHoo < 7 for all a £ A if 
and only if there exists a positive definite polynomial matrix P(a) £ R[a]( n+n ) x ("+" ) such that the inequality in Slip 
is satisfied for all aei. 

Proof: Follows from Theorem 1 in together with Lemma 13.11 above. ■ 



Remark 3.3 To check the condition in Lemma \3.2l first, we should pick an integer dp £ N and search over the set 
of all positive definite polynomial matrices P(a) £ R[oi]*-™ +n ) x ("+" ) such that deg(P(a)) < dp, in order to find a 
feasible solution to the inequality in ill]) for all a £ A. Now, since the degree of P(a) is not known in advance, we 
have to start from an initial value (possibly estimated based on intuition from the physical nature of the problem) and 
keep increasing dp until we reach a feasible solution, which exists if the distance \\G(s; a) — G'(s; a') ||oo is less than 
7. Therefore, we should also start with large values for 7 (to ensure the existence of a feasible solution) and then, 
decrease 7 accordingly (for instance, using the bisection method WN). 



In the next theorem, we use sum-of-squares optimization to rewrite the inequality in (|llj) as a sum-of-square 
feasibility problem which we use later to develop our numerical algorithm. 

Theorem 3.4 Assume that the compact set A can be characterized as 

A = {a £ W I q e (a) > 0, VI < I < L} , 

where qi £ R[a] for all 1 < £ < L. Furthermore, assume that there exist wg £ S[a] for all < £ < L, such that {a £ 
W J wo(a) + ^2i =1 qi(a)wi(a) > 0} is a compact set. Then, for a fixed G'(s;a'), we have \\G(s; a) — G'(s; a')lloo < 7 
for all a £ A if there exist a constant e £ R>o, polynomial matrices P(a) £ iS[a] n+n , and Qi(a) £ <S[a]™ + ™ + m +° 
for all 1 < I < L, such that 



A(a) T P(a) + P(a)A(a) * * 
B(a) T P(a) -I * 

C(a) D{a) -7 2 / 



d + Qo + ^2Qt{a)q e (a) = 0. 

1=1 



(12) 



Proof: The proof follows from Theorem 2 in [16] together with Lemma 13721 above 



Remark 3.5 To check the condition in Theorem \3.7\ we should pick the polynomial degree dp £ N and search over 
the set of all sum-of-square polynomial matrices P(a) £ iS[a] < ' n+n ) x ( n + n ) such that deg(P(a)) < dp, in order to 
solve the polynomial equation in tl3\) . This search is easy to perform since the underlying problem is convex (due to 
the restriction to the set of sum-of-square polynomial matrices) and can be readily solved using available LMI solvers. 
Unfortunately, if we cannot find any solution to this problem for a given degree dp, we cannot deduce that our problem 
does not admit a solution for this given degree dp, since Theorem \S. 7\ is only a sufficiency result. We can only hope 
to find a solution by increasing the polynomial degree dp. 

Note that so far, we assumed that the model matrices A'(a'), B'(a'), C'(a'), and D'(a') are given since otherwise, 
Thcorcm l3.7l would result in nonlinear equations in terms of unknown polynomial coefficients. We propose Procedure[T] 
for finding matrices A'(a'), B'(a'), C'(a'), and D'(a') based on the method of alternating LMIs for solving BMIs [TS]. 

Remark 3.6 This method does not guarantee convergence for the proposed algorithm because there exists always at 
least one BMI which we cannot check its feasibility using the method of alternating LMIs U5f . 

3.2 Complexity Reduction for Discrete-Time Systems 

In the next theorem, we present a result which is a counterpart to Theorem 13.71 for discrete-time systems. 

Theorem 3.7 Assume that the compact set A can be characterized as 

A = {a £ W I q e (a) > 0, VI < I < L} , 

where qe £ R[a] for all 1 < I < L. Furthermore, assume that there exist Wi £ S[a] for all < £ < L, such that {a £ 
R p J wo(a) + 53^=i Qe{ a ) w e( a ) > 0} is a compact set. Then, for a fixed G'(y, a'), we have \\G($; a) — G'(y, a')lloo < 7 
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for all a £ A if there exist a constant e G R>o, polynomial matrices P(a) G S[a] n+n , and Qt{a) G S[a] 2 ( n+n )+ m +° 
for all 1 < I < L, such that 

L 

-eI-Q -Y J Qi{a)qi(a) = 0. (13) 

1=1 

Proof: The proof follows the same reasoning as in Lemmas I3.1H3.2I and Theorem 13.71 ■ 

We can construct a similar procedure for discrete-time systems as in Procedure [T] by changing the nonlinear 
equations in (|P.1[) - (|P.2[) with the nonlinear equation in (fT"3"|) to calculate the reduced discrete-time system. 



P(a) * * * 

P(a)A(a) T P(a) * * 

B(a) T / * 

C(a)P(a) D(a) -f 2 I 



4 Illustrative Example 

In this subsection, we illustrate the applicability of the developed procedure on two numerical examples. The first 
numerical example is a parameter-dependent discrete-time linear systems. We use this example to compare our 
developed algorithm with the method described in |10j . The second example is a parameter-dependent continuous- 
time linear system motivated by controlling power systems. To implement Procedure [TJ we used SOSTOOLS which 
is a free MATLAB® toolbox for formulating and solving sum-of-squares optimizations |25j . 

4.1 Discrete-Time Systems 

Consider the parameter-dependent discrete-time linear system described by 



u(k), 



xi(k + 1)~ 




' 0.5ai 0.1 


x 2 {k + l)_ 




0.3 0.5a 2 




y(k) =[10] 





xi(k) 




' 1 " 




_x 2 (k)_ 


+ 






and 

xi{k) 
x 2 {k) 

where x(k) £ R 2 and u(k) £ R are the state vector and the control input, respectively. Let us define the parameter 
vector as a = [ a\ a 2 ] T G A C R 2 with 

A = {a G R 2 | qi(a) = 1 - a\ > Q 7 q 2 (a) = 1 - a\ > 0} . 

We are interested in reducing the system complexity by getting a new model which is only a function of a±. First, let 
us present the model reduction algorithm introduced in [10] . To do so, we need to introduce the following notations 



where 



and 



Ac 



A(a) = A + axAxUx + a 2 A 2 U 2 , 

, A x = [ 0.5 ], A 2 = [ 0.5 ], 
Ux = [ 1 ] T , U 2 = [ 1 ] T . 



0.1 
0.3 



Now, using [55], it is evident that 

G(r,a) =C{ l I-A(a))- l B 



A, 


Ux 


u 2 


So 


Ax 











A 2 











Co 












3^2X2 

ax 
a 2 



where * denotes the upper linear fractional transformation operator (see |101 .26i for its definition). Let us introduce 
notations 

T 





' A 


Ux 


u 2 ' 




' B Q ' 






A = 


Ax 








, B = 





. c = 







A 2 



















G 



Procedure 2 



Input: Threshold e 6 R>o- 
Output: X,Y G W. 
Initialization: X \d,Y ia G W. 
for k = 1,2,... do 

Solve the optimization problem 

mm x ,Yaw trace(X i d Y + XY \d), 
subject to A T X_A-X + C T _C < 0, 
AY A 1 - Y + BB T < 0, 

if ||X-X old || + ||Y-Y oId || <ethen 

break 
end if 

-^old ^- X. 

Fold <- Y. 
end for 



To get the reduced system, we need to solve the optimization problem 



where 



min j^ygyy trace(JTY), 
subject to AJXA -X + C T _C < 0, 
AYA T - Y + BB T < 0, 



W={W e | W = diag(Wn, W 22 , W 33 ) such that W u E S 2 +1 W 2 2 G 5+ , W 33 e 5^, }. 



(14) 



We use Procedure [2] for solving the optimization problem in (fT4")l . Using [21], we know that if the procedure is 
initialized correctly (i.e., close enough to the optimal solution), the algorithm converges to the optimal decision 
variables. Now, using matrices X, Y G W, we introduce the change of variable T = diag(To, Ti, T 2 ) to get the 
balanced realization of the system 



CT = [ C T ] 



Let us for the moment focus on just removing parameter a 2 from the model matrices (and not decreasing the order 
of the system). Using balanced truncation, we can calculate the reduced system as G r (y, oi\) = C r (}I — A r (ai))~ 1 B r , 
where 





' T ~ 


%T 


Tq 1 UxTx 


T^U 2 T 2 ' 




T 1 B$ 


T -l AT = 


Tf 


1 A 1 T 








, T^B = 







. T 2 


1 A 2 T 














A r { ai ) = T^AoTo + cuT^AtU^t 



0.5a! -1.7 x 1CT 1 
-1.7 x 10" 1 



B r = To -1 B = [ -1-0 0.0 ] T , C r = C T = [ -1.0 0.0 ] 
Finally, we can calculate the error caused by the parameter reduction as 



max ||G r ( 3; ai ) - G( 3 ; a)^ = 0.14 < 2y/a(X 33 Y 33 ) = 0.62, 

where the upper bound of this error was introduced in |10j . 

Now, we can illustrate the result of our proposed algorithm on this numerical example. Let us fix the polynomial 
degrees <1a = 1, = 1, dc = 0, do = 0, dp = 2, dQ = 2, dq 1 = 0, and dQ 2 — 0. We use Procedure Q] when adapted 
for discrete-time systems to get the optimal reduced system with n' = 2. The resulting reduced system is 



5.0 x lO- 1 ^! 
-3.3 x 1CT 1 



A'(ai) = 
For this reduced system, we have 



-1.2 x 1CT 1 
6.5 x 10~ 4 ai 



,B'(ai)=[l.O 6.3 x l(T 2 ai ] T ,C = [ 1.0 ],£»' = 7.9 x 10" 



max ||G'(2;ai) - G(y,a) 



0.095. 
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As we can see, for this particular example, we could achieve a smaller distance between the transfer functions of the 
original system and its reduced one. 

We can also try to reduce the system order by choosing n' = 1. We use Procedure!]] when adapted for discrete-time 
systems to get the optimal reduced system as 

A'{a x ) = 5.2 x lO'V - 2.0 x 1(T 8 , B'{ ai ) = 9.4 x 10"% + 9.9 x 1CT 1 , C = 1.0, D' = 1.8 x 10~ 8 . 

For this reduced system, we have 



max||G'(3;ai) -G(y,a)\\, 

aEA 

while if where using the method in |10j , we would have recovered 



= 0.19, 



max ]|Gy(3; an) - G(y,a)\\oo = 0.27, 

where G r (y, ai) = C r (%I — A r (aij)~ 1 B r with 

A r {ai) = 5.0 x 10 _1 ai, B r ( ai ) = -1.0, C r = -1.0. 



4.2 Continuous-Time Systems 

In this subsection, we present a practical continuous-time numerical example. Let us consider a simple power network 
composed of two generators (see Figure [1]). We have partially extracted the structure of this example and its nominal 
numerical values from |28j . We can model this power network as 

&(t) =Wl(f), 



and 



wi(t) = ^[Pi(t)-iu 1 Bm(<yi(t)-<y a (t))-xr 1 M*i(*))--DiWi(*)] : 

JVL\ 



5 2 (t) = u) 2 (t), 



co 2 (t) = -L [P 2 (t) - x^ 1 sm(S 2 (t) - Sx(t)) - x^ 1 sm(S 2 (t)) - D 2 u 2 {t)] , 

M 2 



where 6i(t) is the phase angel of the terminal voltage of the generator i, u>i(t) is its rotation frequency, and Pi(t) is 
mechanical input power to the generator. We assume that Pi(t) = 1.6 + u\{t) and P 2 (t) = 1.2 + u 2 (t), where u\(t) 
and u 2 (t) are the control inputs to the system. The power network parameters can be found in Table [T] (see |28) for a 
description of these parameters). Note that all these values are given in per unit. Now, we can find the equilibrium 
point (5*, 5|) of these nonlinear coupled systems and linearized the overall system around its equilibrium which would 
result in 
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-X^2 cos C^l — ^2 ) — x l 1 COs((5^) 
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COs(#2 — ) 



Pi 
Mi 








cos(j£--j|) 
x 12 Mi 



s(S2—S^)—X2 1 cos((5^ ) 
W 2 






1 

-D2 

Mi 









"00" 




Aui(t) 




1 




AS 2 (t) 


+ 







Aw 2 (t) 




1 



u 2 {t) 



(15) 



where A8i(t), A5 2 (t), Aui(t), and Aw 2 (t) denote the deviation of the state variables 8i(t), 8 2 (t), uii(t), and u> 2 (t) 
from their equilibrium points. Let us assume that we have connected impedance loads to each generator locally. 
Hence, the parameters x\ and x 2 can vary over time according to the load profiles. Furthermore, assume that each 
generator changes its input mechanical power according these local load variations. Doing so, we would not change 
the equilibrium point (8*,S 2 ). For this setup, we can model the system as a continuous-time parameter-dependent 
linear system described by 



where 



A(a) 



G{s;a) 





-2.7 x lO 1 ^! - 1.5 x 10 2 


7.8 x 10 1 



x(t) = A(a)x(t) + Bu(t), 
y(t) = Cx(t) + Du{t), 



-2.5 x 10- 1 







9.8 x 10 1 


-2.3 x 10 1 a 2 - 1.2 x 10 2 





1 

-2.0 x 10" 



(16) 
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Figure 1: Schematic diagram of the power network in our numerical example. 
Table 1: Nominal value of power system parameters extracted from 1281 . 



Variable 


Nominal Value (p.u.) 


Mi 


2.6 x 1CT 2 


M 2 


3.2 x 1CT 2 


X12 


4.0 x 


x x 


5.0 x lO" 1 




5.0 x 10" 1 




6.4 x 10~ 3 


D 2 


6.4 x 10- 3 



and 



B=[0 10 0] T , C=[l } D = 0, 



with a — [ai Q!2] T - In this formulation, parameter a.i for i = 1, 2. denotes the deviation of the admittance x~ l from 
its nominal value (see Table [1]). Note that here we have chosen the input-output pair to derive a reduced model for 
the network from the perspective of the first generator. One can try to solve this problem for any other given set of 
inputs and outputs. We assume that 



A = {a E M 2 I 0.1 2 - of > for i = 1, 2} . 

Let us fix the polynomial degrees cIa = 1, ds = 0, dc = 0, djj = 0, dp = 3, dQ a = 2, dQ 1 = 2. and d^ 
Procedure [1] to get the optimal reduced system with n' — 4. The resulting reduced system is 



2. We use 



G'(s;ai) 



x(t) = A'(a 1 )x(t)+ B'u(t), 
y (t) =C'x(t)+D'u(t), 



where 



A'(ai) = 



and 



7.4 x 10 _4 qi - 4.2 
-3.0 x lO^i +4.4 x 10" 1 

-6.6 x 10 _3 ai + 1.5 
5.9 x 10" 2 ai - 1.8 x 10~ 2 



5.5 x l(T 4 ai + 2.6 x 1CT 2 



1.7 x 10 _4 ai + 8.8 x 1(T 3 



1.9 x HTV - 6.4 -4.4 x 10" 2 ai - 2.8 x 10 _1 
-2.3 x 10~ 4 ai + 3.8 x 10" 2 2.8 x 10~ 4 ai - 4.2 x 10 _1 



-7.5 x l(T 4 ai + 7.3 x 1CT 2 
-1.4 x 10 _1 ai - 4.4 x 10" 1 
-1.0 x 10~ 4 ai +4.6 x 10~ 2 



-4.7 x 10" 5 ai 



x 10" 1 4.2 x 10" 3 ai + 1.8 x 10 _1 -5.5 x 10" 3 ai 



8.3 x lO" 1 
(17) 



B' 



9.0 x 10~ 2 
3.8 

-4.9 x 10~ 2 
5.6 x 10" 1 



c 



1.1 x 10- 1 
4.3 x 10" 1 
-3.2 x 10" 2 
1.1 x lO" 1 



D' 



-0.10558. 



For this reduced system, we have 



max \\G'(s; a±) — G(s; a) 



= 1.5 x 10" 



Hence, we could effectively remove the model matrices dependencies on the second subsystem parameter a 2 while 
not drastically changing the first subsystem input-output transfer function. 
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5 Conclusions and Future Work 



In this paper, we presented a powerful procedure for approximating parameter-dependent linear systems with less 
complex ones using fewer model parameters or state variables. To do so, we minimized the distance between the 
transfer function of the original system and its reduced version. We presented a suboptimal method for solving 
this minimization problem using sum-of-squares optimization and the method of alternating LMIs for solving BMIs. 
We developed numerical procedures for both continuous-time and discrete-time system contrary to the available 
result which focused mostly on discrete-time systems. Due to relying on sum-of-square optimization, the developed 
procedures would not scale well with the system dimension and the number of parameters. Possible future research 
could focus on developing a better numerical procedure for dealing with BMIs and studying the convergence properties 
of this numerical approach. Furthermore, we could exploit sparsity patterns or symmetry structures to improve the 
scalability of the sum-of-squarc optimization. 
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